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ABSTRACT 


The  problem  of  spontaneous  ionization  (i.  e.  ,  no  externally  applied  elec¬ 
tromagnetic  fields,  nor  hard  radiation)  in  the  reaction  zone  behind  strong 
normal  shock  waves  in  air  has  been  treated  concurrently  with  the  problem  of 
dissociation  and  vibrational  relaxation.  Through  a  comparison  of  specific 
ionization  rates,  one  may  conclude  that  up  to  a  shock  velocity  of  1  0  km/sec 
(about  30  times  the  speed  of  sound  at  room  temperature)  the  predominant 
electron  production  process  would  be  atom -atom  ionizing  collisions.  This 
would  be  followed  in  an  approximately  decreasing  order  of  importance  by 
photoionization,  electron  impact,  atom-molecule  collisions  and  molecule- 
molecule  collisions.  The  charge  exchange  reactions,  while  not  contributing 
directly  to  the  electron  production  process,  were  found  to  have  a  small  but 
noticeable  indirect  effect  on  the  resultant  electron  density  distribution  at 
some  distance  behind  the  shock  due  to  their  continuous  shifting  of  the  rela¬ 
tive  population  between  atomic  and  molecular  ions  (which  recombine  with 
the  electrons  at  different  rates).  The  specific  rate  constants  for  the  atom- 
atom  processes  required  to  interpret  all  existing  experimental  results  ap¬ 
pear  to  be  consistent  with  a  simple  extrapolation  of  the  low  temperature 
rate  constants  according  to  the  crossing -point  model  of  Bates  and  Massey 
for  atom-atom  ionizing  collisions. 


I.  INTRODUCTION 


In  a  complementary  paper,  Lin,  Neal  and  Fyfe^  have  presented  experi¬ 
mental  results  showing  some  typical  electron  density  profiles  behind  normal 
shock  waves  in  air  at  initial  pressure  0.  02  —  p^  ^  0.  Z  mm  Hg.  and  in  the 
shock  velocity  range  4.  5  Ug  ^  1  km/sec  (or  shock  Mach  number  range 
14  ^  Mg  20).  They  have  also  presented  various  measured  quantities 
which  may  be  used  to  characterize  the  ionization  history  behind  the  shock 
front.  The  essential  result  was  that  within  the  range  of  shock  strength 
covered,  the  characteristic  ionization  distance  was  only  about  10  to  40 
times  the  viscosity  mean  free  path  of  the  undisturbed  gas  ahead  of  the  shock, 
and  that  the  electron  density  appeared  to  "overshoot”  the  equilibrium  value 
by  a  considerable  factor  for  some  distance  behind  the  shock  front  before 
relaxing  back  toward  the  final  equilibrium  value. 

The  problem  of  approach  to  equilibrium  ionization  behind  strong  shock 
waves  in  argon  has  been  extensively  treated  by  Petschek  and  Byron,  ^  and 
'also  by  Bond.  ^  The  general  conclusion  was  that  the  main  part  of  the  ioniza¬ 
tion  history  was  dominated  by  the  electron  impact  process,  even  though  the 
starting  mechanism  which  accounted  for  about  10  percent  of  the  total  ioniza¬ 
tion  was  not  clear.  The  local  ionization  rate  at  any  point  within  the  main 
ionization  zone  was  then  completely  determined  by  the  local  degree  of  ioni¬ 
zation  (i.  e.  ,  fraction  of  atoms  ionized)  and  the  electron  temperature.  The 
electron  temperatufe  was  determined  by  a  balance  between  the  rate  at  which 
electrons  lose  energy  in  the  ionizing  collisions  and  the  rate  at  which  they 
regain  energy  from  elastic  collisions  with  the  hotter  positive  ions  and  the 
neutral  atoms.  The  translational  temperature  of  the  latter  was,  in  turn, 
uniquely  determined  by  the  local  degree  of  ionization  and  the  energy  equa¬ 
tion.  ^  The  resultant  effect  was  that  the  electron  density  increased  mono- 
tonically  in  an  approximately  exponential  manner  until  the  backward  rate 
(i.  e.  ,  electron-ion  recombination)  became  sufficiently  rapid  to  limit  the 
net  rate  of  electron  production. 

While  the  approach  to  ionization  equilibrium  behind  extremely  strong 
shock  waves  in  air  (i.  e.  ,  kinetic  energy  per  air  molecule  much  greater 
than  the  dissociation  energy)  may  assume  some  of  the  characteristics  of 


the  argon  shock  just  described,  the  predominant  ionization  path  is  likely  to 
be  very  different  at  more  moderate  shock  velocities  where  the  dissociation 
energy  remains  an  important  part  of  the  total  energy  storage  behind  the 
shock.  The  presence  of  more  than  one  chemical  species  in  the  initial  gas 
mixture  is  also  expected  to  further  complicate  the  air  problem. 

Results  obtained  from  studies  of  relaxation  phenomena  in  air  and  its  com¬ 
ponent  gases  have  been  reported  by  the  authors  and  their  colleagues  on  vari¬ 
ous  occasions  over  the  past  years.  These  studies  have  included  critical 

examination  of  the  steadily  growing  volume  of  chemical  kinetic  data,  but  the 
ionization  processes  have  only  been  treated  in  a  cursory  manner. 

In  the  present  paper,  the  various  ionization  processes  that  may  play  sig¬ 
nificant  roles  behind  strong  shock  waves  in  air  will  be  critically  examined. 
The  composite  ionization  history  behind  the  shock  will  be  constructed  from 
these  processes  along  with  the  chemical  and  vibrational  relaxation  processes 
for  various  shock  conditions.  The  theoretical  results  will  then  be  used  to 
interpret  the  experimental  results  reported  in  the  complementary  paper. 
Some  remarks  will  also  be  made  concerning  the  validity  of  extrapolating  the 
present  theoretical  model  to  a  wider  range  of  shock  conditions. 

II.  COUPLING  BETWEEN  VIBRATIONAL  EXCITATION, 
DISSOCIATION  AND  IONIZATION 

It  is  well  known  that  the  combined  problem  of  vibrational  and  radiative 
excitations,  dissociation,  and  ionization  behind  strong  shock  waves  in  a  dia¬ 
tomic  gas  mixture,  such  as  air,  is  very  complex  on  account  of  the  large 
number  of  species  and  reactions  which  must  be  simultaneously  considered 
and  followed.  ^  However,  in  a  situation  where  the  degree  of  ionization 
remains  sufficiently  low  and  the  radiative  heat  loss  is  small  compared  with 
the  rate  of  enthalpy  flow  across  the  shock  front,  one  may  take  advantage  of 
the  "one-way  coupling"  between  the  dissociation  process  and  the  ionization 
and  radiative  excitation  processes.  In  other  words,  while  the  local  ioniza¬ 
tion  and  radiative  excitation  rates  at  any  given  point  behind  the  shock  front 
are  sensitive  to  the  local  temperature  and  chemical  composition  resulting 
from,  the  vibrational  excitation  and  molecular  reactions,  they  do  not  con¬ 
tribute  significantly  to  the  reaction  rates  of  these  molecular  processes, 
nor  do  they  affect  the  energy  balance  in  a  significant  manner.  In  such  a 
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situation,  one  may  first  determine  the  temperature,  density,  and  chemical 
composition  profiles  independent  of  the  ionization  and  photo  processes.  From 
these,  one  can  then  calculate  the  ionization  and  radiation  profiles  from  the 
corresponding  specific  rates  for  these  processes.  However,  at  higher 
shock  strengths,  such  an  iterative  procedure  would  not  be  satisfactory,  and 
one  must  consider  the  electronic  and  photo  processes  simultaneously  with 
the  chemical  processes  in  the  shock  profile  determination.  For  the  case  of 
air,  the  dividing  point  for  the  two  situations  should  be  somewhere  about  10 
km/s  ec  shock  velocity,  where  the  maximum  degree  of  ionization  behind  the 
shock  does  not  appreciably  exceed  one  percent. 

In  the  present  paper,  we  shall  consider  only  the  simpler  situation  in  which 
the  coupling  between  dissociation  and  ionization  is  approximately  one  way. 
Then,  considering  air  to  be  an  initially  binary  mixture  consisting  of  21  yo 
and  79%  N2  by  partial  pressure  (inclusion  of  argon,  CO^,  and  other 

minor  constituents  would  only  complicate  the  computational  work  without 
contributing  to  the  essential  results),  one  may  proceed  to  determine  the 
temperature,  density,  molecular  and  atomic  concentration  profiles  behind 
the  shock  from  the  following  major  chemical  reactions; 
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The  chemical  symbols  employed  in  the  above  expressions  are  conventional, 
except  that  we  have  used  X  to  represent  any  molecule  or  atom  in  the  disso¬ 
ciating  gas  mixture.  The  energy  of  reaction  at  0°K  is  given  in  electron- 
volts  (e.  V.  ).  Reactions  (1)  through  (5)  are  the  same  as  those  considered  in 
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earlier  works  by  the  authors^’ 8,10  others^  in  connection  with  the 

problems  of  chemical  reactions  behind  shock  waves  and  in  the  flow  fields  of 
high  temperature  air.  The  bimolecular  reaction  (6)  was  thought  to  be  unim¬ 
portant,  but  recent  studies  of  the  decomposition  rate  of  NO  behind  shock 
waves  by  Freedman  and  Daiber  and  by  Wray  and  Teare  suggested  that 
the  rate  constant  for  this  reaction  in  the  temperature  range  of  interest  could 
be  much  greater  than  what  had  been  previously  estimated.  ^  ^ This  reac- 

tion  has  subsequently  been  included  in  more  recent  works.  ’  It  was 
found  that  this  reaction  does  affect  the  initial  rate  of  formation  of  NO  behind 
the  shock,  but  the  overall  effect  on  the  ionization  processes  appears  to  be 
minor. 

From  the  set  of  chemical  reactions  shown  above,  and  from  the  usual 
definition  of  the  constant-volume  reaction  rate  constant  K ,  k  (super¬ 
script  referring  to  the  catalyst  X),  etc.  ,  one  may  write  for  the  total  time 
rate  of  change  of  molecular  oxygen  number  density  following  any  fluid  element. 


ri  02*^  X  ■  '^Rl  ^>0*^ 

'^5'^N^02  ^  '^RS'^NO'^O  (7) 

'^6'^N2'^02  ■'■'^Rb'^NO 

Or,  in  terms  of  the  normalized  number  density  [,X]Env/n  •  ,  where 
n  .  =  (p/p,  )  n  ,  is  the  equivalent  number  of  air  molecules  per  unit 

volume  at  the  local  mass  density  p  behind  the  shock,  and  H  ,  Pj^  denote, 
respectively,  the  initial  number  density  and  mass  density  ahead  of  the 
shock. 
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IX]  +  Kj  [NJ  +  [N3] )  [O^J  n  1 

X  1 

+  (Kj^j  (NO]  [O]  +  [NO]^)  ^  n  J  (7-a) 

+  [o]^[x]  rii^ 


Similar  equations  can  be  written  for  the  total  time  rate  of  change  of  the 
number  density  for  other  chemical  species. 

In  addition  to  the  chemical  reactions  listed  above,  one  must  also  con¬ 
sider  the  effects  of  molecular  vibrational  relaxation,  and  of  the  coupling 
between  vibrational  relaxation  and  dissociation.  As  has  been  proposed  by 
Hammerling,  Teare  and  Kivel,"^  these  effects  can  be  approximately  taken 
into  account  by  applying  a  correction  factor  of  the  form 


1  -  exp  [  -N  (0  -  ©)  1 

_  V  V 

[exp  (0^  -  0)  -  1] 


exp  0  -  1 

_ V 

exp  0-1 


(8) 


to  the  dissociation  rate  constant  corresponding  to  local  vibration-translation 
equilibrium  (k  .  such  that 

■  V  C?) 

In  Eq.  (8),  denotes  the  total  number  of  vibrational  levels  for  the  molecule 
considered  (assuming  to  be  finite  and  the  levels  equally  spaced);  0^  =  Ev/ 
kT^  and  0  =  E.^^/kT  are  respectively  the  reciprocal  of  the  vibrational  and 
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translational  temperature  normalized  with  respect  to  the  vibrational  energy 
spacing  (k  being  the  Boltzmann  constant).  The  vibrational  temperature 
is  related  to  the  translational  temperature  through  the  usual  relaxation  ex¬ 
pression 


d£ 

V 

"Ht 


(£  ) 

V  eg. 


(10) 


where  =  E^/(exp  0  ^  -  1)  is  the  mean  vibrational  energy  per  molecule; 
and  is  the  vibrational  relaxation  time,  which  may  be  assumed  to  depend 
only  on  the  translational  temperature  and  particle  density.  18,  19 

The  coupling  mechanism  between  vibrational  relaxation  and  dissociation 
discussed  above  is  equivalent  to  assigning  equal  weight  to  the  pre -exponential 
factor  for  the  partial  dissociation  rate  constant  for  the  various  vibrational 
states,  and  it  assumes  that  the  population  of  the  vibrational  states  fol¬ 
lows  the  Boltzmann  distribution  at  all  times  throughout  the  relaxation 
process. 

The  conservation  equations  for  a  plane  shock  wave  propagating  at  a  con¬ 
stant  velocity  Ug  in  an  otherwise  undisturbed  atmosphere  of  density  and 

pressure  p,  are  given  by  (neglecting  viscous  dissipation  within  the  reaction 

^  21 
zone  and  radiative  heat  loss). 


Mass; 


pu  =  p  U 
1  s 


(11) 


Momentum: 


p  +  pu 
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+  P, 
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Energy: 
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P, 
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2 

u 

2“ 


£l 

m 


(13) 


where  u  is  the  local  mass  velocity  in  a  coordinate  system  fixed  to  the  propa¬ 
gating  shock  front,  p  =  kT  S  D  is  the  local  gas  pressure,  and 

X  X 


f  =  E  [X]  f  (14) 

X  X 


is  the  total  thermal  and  chemical  energy  per  equivalent  air  molecule,  and 
m  is  the  mass  per  air  molecule.  Note  that  C  depends  only  on  the  local 
translational  temperature  T,  the  vibrational  temperature s  (T^)02,  (Tv)n2’ 
and  the  normalized  particle  densities  |]x]  . 

By  integrating  the  chemical  rate  equations  (of  which  Eq.  (7)  is  a  typical 
example)  and  the  vibrational  relaxation  equations  for  O2  and  N2  simultaneously 
under  the  constraints  of  these  conservation  equations,  one  can  follow  the  trans¬ 
lational  and  vibrational  temperatures,  mass  density,  and  particle  densities 
as  a  function  bf.time  t,  or  distance 

.3r  =  /*  udt 

behind  the  translational/rotational  shock  front  for  any  given  set  of  initial 
conditions.  This  is  illustrated  in  Figs.  1(a),  (b),  (c)  and  2  (a),  (b),  (c)  where 
the  local  temperatures,  mass  density,  and  chemical  composition  behind  a 
plane  shock  wave  in  air  at  three  different  velocities  in  the  range  of  initial 
pressure  0.  02  <  2  mm  Hg  are  plotted  as  functions  of  the  normalized 

distance  behind  the  shock  (i.  e.  ,  normalized  with  respect  to  the  upstream 
mean  free  path  j2  j^).  The  vibrational  temperature  for  NO  can  be  assumed 
to  be  identical  to  the  translational  temperature  since  this  chemical  species 
is  produced  only  through  energetic  collisions  behind  the  shock. 

The  method  of  solution  and  computational  procedure  employed  in  obtaining 
the  numerical  results  presented  here  are  essentially  the  same  as  those  des¬ 
cribed  in  Ref,  12.  It  may  be  mentioned  again  that,  as  a  simplifying  approxi¬ 
mation,  vibrational  relaxation  and  chemical  reactions  within  the  translational/ 
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rotational  relaxation  zone  have  been  completely  neglected,  so  that  the  start¬ 
ing  point  X  =  O  actually  refers  to  some  fictitious  point  at  which  the  gas  tem¬ 
perature  and  density  assume  values  corresponding  to  full  translational/rota- 
tional  relaxation  behind  the  shock.  Since  the  actual  characteristic  distance 

for  translational/  'otational  relaxation  behind  strong  shock  waves  is  known 

6  Z  Z 

to  be  comparable  to  the  upstream  mean  free  path,  '  the  physical  meaning 
for  the  region  X  / Jl-^  ^1  in  the  present  model  becomes  ambiguous.  For  the 
region  X  /-^i  >  1,  there  would  also  be  an  uncertainty  of  the  order  of  unity  in 
the  absolute  scale  of  X  with  respect  to  the  center  of  the  translational/ 

rotational  relaxation  zone.  These,  however,  may  no^  be  considered  as  seri¬ 
ous  limitations  as  long  as  the  characteristic  distances  for  vibrational  relaxa¬ 
tion  and  chemical  reactions  are  still  considerably  greater  than  ^ 

The  chemical  rate  constants  employed  in  the  present  calculations  are 
summarized  in  Table  I.  Aside  from  the  addition  of  the  bimolecular  reaction 

(6),  this  set  of  rate  constants  is  essentially  the  same  as  those  quoted  in 

N-j 

Refs.  4  and  6,  with  the  exception  of  an  upward  revision  of  K  and 

K  L-j  ,  by  a  factor  of  2  and  4,  respectively;  a  reassignment  of  the  cata- 
lytic  efficienty  of  O  and  N  in  reaction  (3),  and  some  minor  adjustments  of 
the  temperature  dependence  for  and  according  to  the  recent  re¬ 

view  paper  by  Wray.  '  The  equilibrium  constants  employed  to  relate  the 
backward  and  forward  rate  constants  for  each  reaction  are  based  on  standard 
spectroscopic  data.  23-25 

The  most  striking  features  of  the  temperature  and  chemical  composition 
profiles  for  strong  shock  waves  in  air  as  exhibited  in  Figs.  1  and  2  are 
perhaps  the  rapid  progress  of  the  oxygen- nitrogen  reactions  immediately 
behind  the  shock  front  on  the  one  hand  and  the  relatively  long  time  it  takes 
some  of  the  reactions  to  approach  the  final  adiabatic  equilibrium  state  be¬ 
hind  the  shock  on  the  other  hand.  For  example,  at  a  distance  of  10  mean 
free  paths  behind  the  translational/rotational  shock  front  for  the  median  shock 

*It  may  be  mentioned  that  Wray  in  Ref.  15  has  suggested  adopting  the  same 
temperature  dependence  for  all  catalyst  X  in  reaction  (1).  However,  this 
would  tend  to  overestimate  the  relative  catalytic  efficiencies  for  O  and  O2 
at  high  temperatures  (i.  e.  ,  T  >  8000  K).  Therefore,  the  rate  constants 

/C  and  as  quoted  in  Refs.  4  and  6  are  still  preferred. 
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TABLE  I 


Summary  of  chemical  rate  constants  employed  in  neutral  particle 
concentration  determination  (T  in  °K,  concentration  in  particles/ 

3 

cm.  ,  time  in  seconds). 


RATE  CONSTANT 

CATALYST  X 

1  6.  2  X  10"^® 

O 

R1 

I  2.  2  X  10"^® 

^2 

1  1.  7  X  10"^^ 

^2 

‘  8.  3  X  10"^^ 

N,  NO 

1 

6.  5  X  lO"^"^ 

N 

kX  =  1 

7.  6  X  10-^2  T-'/^ 

^2 

R2  1 

3.  0  X  10'^^ 

O2,  0,  NO 

K  X  = 

R3 

5.  5  X  lO"^”^ 

NO 

2.  8  X  10"^® 

0,  N,  O2,  N^ 

'^R4  = 

2.  7  X  10'^^ 

'^5  " 

2.  2  X  10"^^  T  exp  (-3,  56o/T) 

R6 

0.  4  T'^/^  exp  (-43,  OOO/T) 
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velocity  cited  (i.  e.  ,  =  6.  9  km/sec),  dissociation  of  and  have 

reached  7  5%  and  50%  of  the  final  equilibrium  level  respectively,  and  the 
translational  temperature  has  dropped  by  a  factor  of  two  from  the  initial 
value.  On  the  other  hand,  it  takes  a  distance  of  several  hundred  mean 
free  paths  (depending  on  the  initial  pressure  pj^ )  in  order  for  dissocia¬ 
tion  and  the  translational  temperature  to  come  within  10%  of  the  final  equi¬ 
librium  level.  The  rapid  initial  rates  of  change  in  temperature  and  chemi¬ 
cal  composition  imply  that  the  various  competing  ionization  processes  may 
carry  different  weights  at  different  distances  behind  the  shock  front.  The 
long  chemical  relaxation  time  toward  the  equilibrium  state  implies  that  the 
degree  of  ionization  may  also  be  expected  to  depart  from  the  final  equilib¬ 
rium  level  for  a  large  distance  behind  the  shock  front. 

Further  discussion  of  the  vibrational  and  chemical  relaxation  histories 
is  beyond  the  scope  of  the  present  paper,  but  an  illustration  of  how  the  vari¬ 
ous  reactions  proceed  behind  the  shock  front  can  be  found  in  the  recent  paper 
1 5 

by  Wray.  Similar  illustration  can  also  be  found  in  an  earlier  paper  by  Duff 
g 

and  Davidson,  in  which  substantially  different  rate  constants  were  used  for 
some  of  the  important  reactions.  The  effect  of  this  different  choice  of  the 
rate  constants  were  discussed  in  Ref.  6. 

III.  THE  COMPETING  IONIZATION  PROCESSES 

To  determine  the  ionization  history  behind  the  shock  front  in  the  energy 
range  of  interest  (see  preceding  section)  we  shall  consider  the  following 


classes  of  ionization  processes: 

Electron-Impact  and  Ion-Impact  Ionization 

e  +  X  Ijrzr  X+  +  e  +  e  (15) 

X'''  +  Y  X'*'  +¥''■+  e  (16) 

Photoionization 

hv  +  X  <  **  X'*'  +  e  (17) 

Charge  Exchange 

X'''  +  Y  X  +  Y'''  (18-a) 

A'^  +  ^  AB  +  B'^  (18-b) 

X'*'  +  Y  ^  -%  X  +  Y"*”^  +  e  (18-c) 
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Electron  Attachment 


e 

+ 

■5^  ■; — 

:x 

(19-a) 

e 

+ 

~ — 

:  X  +  h.v 

(19-b) 

e 

+ 

AB  - - ► 

:A  +  B“ 

(19-c) 

Ionization  by  Neutr 

al  Atom  and  Molecule  Collisions 

X 

+ 

Y  - - ► 

X'*'  +  Y  +  e 

(20) 

A 

+ 

R  - ► 

AB'*'  +  e 

(21) 

A 

+ 

B2=^ 

AB^"*"  +  e 

(22-a) 

A 

+ 

B2^ 

AB'*'  +  B  +  e 

(22-b) 

^2 

+ 

^2  — 

A^B  +  B  +  e 

(23-a) 

^2 

+ 

AB'*'  +  AB  +  e 

(23-b) 

In  the  above 

expr e  ssions. 

we  have  used  the  symbols,  X,  Y 

,  to  represent 

any  neutral  molecule  or  atom;  the  symbols  A,B,  to  represent 

a  neutral  atom 

only  (,i.  e.  ,  diatomic  molecules  to  be  represented  by  A^,  AB,  etc.  ),  and 

the  symbol  e  to  represent  a  free  electron.  The  above  list  is  selective  rather 
than  comprehensive,  in  the  sense  that  only  processes  involving  the  lowest 
energies  of  reaction  within  each  group  have  been  included.  We  have  also 
excluded  electron  production  processes  requiring  3-body  collisions,  since 
they  would  be  relatively  infrequent  in  the  range  of  gas  density  under  con¬ 
sideration.  To  avoid  further  complication  by  the  question  of  distribution 
among  the  various  vibrational  and  electronic  states  for  the  molecules,  atoms, 
and  ions  involved  in  Eqs.  (15)  through  (23-b),  we  shall  henceforth  regard  all 
the  cross-sections  (or  rate  constants)  associated  with  these  electronic  proc¬ 
esses  as  averaged  quantities  over  the  appropriate  internal  energy  states,  un¬ 
less  otherwise  noted. 

In  addition  to  the  gas-phase  reactions  listed  above,  one  may  also  consider 
various  secondary  processes,  such  as  electron  diffusion  and  photo-electric 
effects  from  solid  surfaces,  which  have  been  observed  in  some  noble  gas 
experiments.  26,27  Since  these  effects  are  generally  weak  compared  v/ith 
the  main  ionization  processes,  and  are  to  some  extent  geometry-dependent, 
we  shall  postpone  our  discussion  of  these  effects  until  later  in  Section  IV. 
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Let  us  consider  now  the  various  gas -phase  ionization  processes  in  some 
detail.  In  view  of  the  fact  that  the  forward  and  backward  rates  are  closely 
related  for  each  process  through  the  local  equilibrium  constant,  we  only 
need  to  discuss  either  one  of  the  two  rates  in  each  case. 

A.  Electron-Impact  and  Ion-Impact  Ionization 

Electron-impact  ionization  has  long  been  known  to  be  the  single  most 
important  electron  production  process  in  almost  every  gas  discharge  phe¬ 
nomenon.  It  has  been  found  to  be  also  the  main  ionization  process 

2  3 

behind  strong  shock  waves  in  noble  gases.  > 

In  terms  of  the  local  electron  energy  distribution  function  f  (e),  and  the 
total  ionization  cross  section  collisions  between  electrons  of 

energy  e  and  molecules  or  atoms  of  type  X,  the  instantaneous  rate  of  pro¬ 
duction  of  ion  pairs  per  unit  volume  following  a  fluid  element  behind  the 

30 

shock  front  is  given  by  (considering  only  the  forward  rate). 

V  (e)  .  f  (£)  d  e  (24) 

where  n  is  the  local  number  density  of  free  electron;  w^x  threshold 

energy  for  the  ionization  cross  section  Q^X’  v  =  (2  e/m^)^/  is  the  electron 
speed  corresponding  to  energy  e  (m^  being  the  electron  mass). 

In  the  absence  of  strong  electromagnetic  fields,  the  electron  energy  dis¬ 
tribution  will  be  approximately  Maxwellian  at  a  temperature  Tg  which  is 
comparable,  but  not  necessarily  equal,  to  the  local  translational  gas  tem¬ 
perature  T.  Strictly  speaking,  the  electron  temperature  profile  behind  the 
shock  front  should  be  determined  from  a  detailed  consideration  of  the  local 
energy  balance  for  the  free  electrons  within  each  fluid  element,  taking  into 
account  the  effect  of  heat  conduction  by  the  electrons  (due  to  the  high  thermal 
conductivity  of  the  electrons  and  the  large  variation  in  temperature  within 
the  reaction  zone);  as  well  as  the  effects  of  energy  transfer  between  the 
heavier  particles  and  the  electrons;  energy  loss  by  the  electrons  in  impact 
ionization  and  excitation,  etc.  ^  However,  in  view  of  the  fact  that  energy 
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transfer  for  collisions  between  diatomic  molecules  and  electrons  appears 
to  be  quite,  efficient,  *  and  that  the  electron  mean  free  path  is  generally 
much  smaller  than  the  thickness  of  the  reaction  zone,  the  actual  lag  between 
the  two  temperatures  may  not  be  too  great.  Therefore,  in  a  first  approxi¬ 
mation,  we  may  consider  Tg  to  be  everywhere  the  same  as  T  in  estimating 
the  impact  ionization  rate.  Thus,  Eq.  (24)  becomes 


(15) 


s  n 


7rm 


1/2 

-)  (• 


3/2 


kT 


£ 

X 


(£)  ede 
(25) 


for  kT  <<  the  above  equation  may  be  approximated  by  its  asymptotic 

equivalent, 


'^iX  ^'iX  ('^iX^ 


(26) 


where  Q’  is  the  slope  of  the  Qjx  ^  curve  at  (or  near)  the  thres¬ 

hold  energy 

The  ionization  efficiencies  by  electron  impact  for  N2,  O^,  and  NO  and  a 
number  of  other  molecular  gases  at  room  temperature  have  been  measured 
by  Tate  and  Smith.  32  The  apparent  ionization  cross-section  as  deduced  by 
Massey  and  Burhop**  from  these  measurements  is  reproduced  here  in  Fig. 

3  for  electron  energies  up  to  several  times  the  threshold  energy.  It  is  seen 
that  these  cross-sections  are  all  characterized  by  an  approximately  linear 
rise  near  the  threshold,  which  continues  on  for  a  considerable  distance 

*See  for  example,  Massey  and  Burhop,  loc.  cit.  ,  Table  VI,  p.  279. 

S.  W.  Massey  and  E.  H.  S.  Burhop,  loc.  cit..  Table  III,  p.  265. 
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beyond  the  threshold;  a  maximum  value  which  is  several  times  the  atomic 
unit  of  area  tt  (where  =  5.  29  x  10  ^  cm  is  the  first  Bohr  radius)  at 

about  100  e,  v.  ;  and  then  a  gradual  fall-off  at  higher  electron  energies.  It 
may  also  be  noted  that  the  products  iX  these  molecules  are 

all  very  nearly  equal  to  . 

For  atomic  oxygen  O,  and  atomic  nitrogen  N,  which  become  two  of  the 
most  abundant  species  beginning  at  some  distance  behind  strong  shock  waves 
in  3.ir,  the  electron  impact  ionization  cross-sections  are  not  very  well  known. 
However,  judging  from  the  fact  that  the  cross-sections  vs  electron  en¬ 
ergy  curves  for  all  the  observed  molecules  are  similar  except  for  a  dis¬ 
placement  in  the  threshold  energy,  and  that  they  are  also  very  similar  to 
those  observed  for  noble  gas  atoms,  *  it  is  reasonable  to  assume  that  the 
electron  impact  ionization  cross-sections  for  O  and  N  atoms  are  also  quite 
similar  to  those  for  N2,  02»  and  NO,  except  for  the  difference  in  threshold 
energy.  Therefore,  in  calculating  the  total  electron  impact  ionization  rate 

using  the  asymptotic  formula  (26),  one  may  simply  take  w.^Q'..^,.  (  W'-y-)  = 

2  iX  iX  lA.' 

TT  a  ^  for  all  the  atoms  as  well  as  for  the  molecules. 

The  adoption  of  Eqs.  (25)  or  (26)  together  with  the  cross-sections  obtained 
from  room  temperature  measurements  may  lead  to  an  underestimate  of  the 
total  electron  impact  ionization  rate  behind  the  shock  due  to  the  implicit 
neglect  of  multistep  ionization  via  the  upper  electronic  states.  ^  However, 
such  underestimate  may  have  already  been  compensated  for  by  our  earlier 
assumption  of  identical  electron  and  gas  temperatures. 

For  ion-impact  ionization  [Reaction  (I6)]  ,  the  cross-sections  are 
generally  much  smaller  than  those  for  electron-impact  ionization  at  energies 
below  a  few  hundred  electron  volts.  **  Considering  the  fact  that  the  heavier 
mass  also  tends  to  discourage  the  frequency  of  ion-neutral  collisions,  it 
becomes  quite  clear  that  the  ion-impact  processes  will  be  insignificant  in 
comparison  with  the  electron-impact  processes  for  producing  ionization 
behind  shock  waves. 


*H.  S.  W.  Massey  and  E.  H.  S.  Burhop,  loc.  cit.  ,  Section  6.  4,  p.  530. 

**H.  S.  W.  Massey  and  E.  H.  S.  Burhop,  loc.  cit.  ,  Chapter  VIII. 
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B.  Photoionization 


The  rate  of  production  of  ion  pairs  per  unit  volume  due  to  photoionization 
following  a  fluid  element  can  likewise  be  expressed  as 

(^)  =  ^  “fx  fin  (-)  d-  (^7) 

'^iX/h 

where  C  is  the  velocity  of  light;  h  is  Planck's  constant  (f)  is  the  photo¬ 

ionization  cross-section  for  type  X  molecule  at  light  frequency  v;  and  n  (v) 
dv  is  the  local  number  density  of  photons  in  the  frequency  interval  between 
V  and  V  +  di*. 

The  far  ultraviolet  absorption  spectra  for  N^,  and  NO  at  room  tem¬ 
perature  have  been  studied  by  Weissler,  Po  Lee,  and  others.  The  photo - 
ionization  cross-section  as  deduced  from  these  spectra^'^  are  reproduced 
here  in  Fig.  4  as  functions  of  the  photon  energy  hv.  It  is  seen  that  the  cross- 
section  generally  rises  steeply  from  the  threshold  to  a  peak  value  of  about 
2  X  10"1^  cm^,  and  then  decreases  more  gradually  at  higher  photon  energies. 
While  the  photoionization  cross-sections  for  the  atomic  species  N,  O  are  not 
known,  they  are  not  likely  to  be  very  different  from  those  of  the  molecular 
species,  (after  allowing  for  the  differences  in  threshold  energies)  as  one  can 
judge  from  the  observed  cross-sections  for  other  monatomic  gases  and 
metallic  vapors. 
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The  peak  value  of  about  2x10  cm  for  the  photoionization  cross- 
sections  implies  that  ionizing  photons  of  energy  substantially  above  the 
threshold  may  have  a  mean  free  path  within  the  reaction  zone  behind  the 
shock  that  is  comparable  to  the  thickness  of  the  reaction  zone  itself.  * 


'i'The  absorption  mean  free  path  for  the  ionizing  photons  at  a  local  gas  den¬ 
sity  p  is  then  given  by  Jl  =  300j2jPj/p  where  and  Pj^  are  respectively 
the  viscosity  mean  free  path  and  the  gas  density  ahead  of  the  shock.  Thus, 
in  the  region  behind  the  shock  front  where  p/pj^  =  10,  we  have  ^  30ij^, 
which  is  comparable  to  the  thickness  of  the  reaction  zone. 
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Thus,  the  local  photon  density  n  (v)  &v  in  the  far  ultraviolet  can  only  be 
determined  from  a  detailed  consideration  of  the  excitation,  emission,  and 
radiative  transport  problems  within  the  reaction  zone.  These  problems 
are  generally  quite  complicated  and  the  basic  information  required  to  treat 
these  problems  is  at  present  only  very  coarsely  known.  We  shall, 

therefore,  only  attempt  to  make  a  rough  estimate  of  the  photoionization  rate 
here. 

F rom  a  study  of  the  photoelectric  measurements  by  Camm,  et  al, 
Hammerling^^  has  tentatively  concluded  that  most  of  the  far  ultraviolet 
radiation  from  the  reaction  zone  behind  strong  shock  waves  in  air  (up  to 
Ug  =  7  km/ sec)  may  be  attributed  to  the  b*  ^2+^  — ►  X  transition  of  N2. 

Neglecting  absorption  and  re-emission,  the  total  intensity  of  this  radiation 
per  unit  area  of  the  shock  front  may  be  expressed  in  the  form 


exp  (  -  E'/kT") 


I  S 


E  L 


2  (t  +  t  ) 
r  c' 


(28) 


where  nj^2  i®  the  averaged  number  density  of  N2  molecules  within  the  reac¬ 
tion  zone;  E*!*  =  12.  9  e.  v.  is  the  energy  difference  between  the  two  electronic 
states  of  N2  (measured  from  their  respective  ground  vibrational  levels);  L 
is  the  effective  thickness  of  the  reaction  zone;  and  are  respectively 
the  radiative  life  time  and  collisional  de-excitation  time  for  the  upper  elec- 
tronic  state  of  N2;  and  T  is  a  temperature  which  characterizes  the  popula¬ 
tion  of  the  electronic  states  of  N2  within  the  reaction  zone  when  t^/t^ — >-0. 

By  letting  Q£^(t')  =  2x10  cm^  =  tt ao/4  for  all  photon  energies  above  the 
ionization  threshold  (i.  e.  ,  hf  >  w^^)  and  by  equating  I/e'*'  to  the  total 
ionizing  photon  flux  density,  one  obtains,  from  Eqs.  (27)  and  (28) 


8  (t  +  t  ) 
r  c 


exp  (  -EVkT*) 


(29) 
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It  may  be  noted  that  in  arriving  at  the  above  expression,  we  have  replaced 

Eq.  (27)  by  (no2  since  photons  of  12.  9  e.  v.  energy  would  only 

be  sufficiently  energetic  to  ionize  and  NO  molecules,  but  not  other  chem¬ 
ical  species  behind  the  shock  front.  For  numerical  estimate  of  this  photo¬ 
ionization  rate,  one  may  take,  following  the  suggestion  of  Hammerling, 

Tj.  =  4x10"9  . 


T 

C 


)■ 


where 

is  the  mean  thermal  speed  of  N2  and  ^  ^  10”^^  cm^  is  the  gas  kinetic 

(or  viscosity)  cross-section  of  N2;  and  L,  to  be  twice  the  distance  behind  the 
shock  front  at  which  T  =  T  .  The  effective  temperature  T',  which  is  by  far 
the  most  sensitive  parameter  in  Eqs.  (28)  and  (29),  is  somewhat  uncertain. 
For  fitting  the  limited  experimental  results  by  Camm,  et  al,  one  may 

take  T  _  (T^  — » o  ^x— >c3o  )  where  ^  is  the  gas  temperature  imme¬ 
diately  behind  the  translational/rotational  shock  front,  and  T  is  the 

X  —►00 

final  equilibrium  temperature  far  behind  the  shock. 

C.  Charge  Exchange 

The  charge  exchange  reactions  (18-a)  and  (18-b)  involve  a  transfer  of 
unit  charge  from  one  molecule  (or  atom)  to  another  during  a  collision, 
and  hence  do  not  contribute  to  the  net  electron  production  rate  during  the 
early  part  of  the  ionization  history.  However,  at  later  times  when  the  elec¬ 
tron-ion  recombination  rate  becomes  comparable  to  the  ionization  rate, 
such  charge  exchange  reactions  could  influence  the  electron  density  profile 
through  a  continuous  reshuffling  of  the  relative  populations  among  the  differ¬ 
ent  positive  ions.  This  effect  would  especially  be  significant  if  there  exists 
a  large  difference  in  recombination  rates  for  the  different  kinds  of  positive 
ions  (e.  g.  ,  radiative  recombination  for  atomic  ions  versus  dissociative 
recombination  for  molecular  ions). 

For  the  oxygen-nitrogen  system  under  consideration,  the  following  speci¬ 
fic  reactions  are  of  the  straight  forward  type  (18-a), 
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N2‘''  +  N  -5— *• 

n/  +  O 

n/  +  O2  5= 

N2‘''  +  NO  :^=r*-  N2 

+  o  -r~*"  N 


+ 

N+ 

+ 

1. 0 

e.  V, 

+ 

+ 

0 

+ 

2.  0 

e.  V. 

+ 

02^ 

+ 

3.  5 

e.  V. 

+ 

no'*' 

+ 

6,  3 

e.  V. 

+ 

+ 

0.  9 

e.  V. 

N"^  +  O2  N 

N"^  +  NO:5=i  N 


+ 

0 

tv 

+ 

+ 

2.  5  e, v, 

+ 

no'*' 

+ 

5.  3  e. V 

o 


+ 


o 


+ 


+  ^2^  1 . 6  e.  V. 

+  no"*"  +  4.  4  e.  V. 


NO: 


O. 


+  no"''  + 


2.  8  e.  V. 


(30) 

(31) 

(32) 

(33) 

(34) 

(35) 

(36) 

(37) 

(38) 

(39) 


For  these  charge  exchange  reactions,  quantitative  information  on  the  rate 

constant  (or  cross-section)  is  very  meager,  especially  at  the  relatively  low 

particle  energy  under  consideration,  30  Existing  data,  however,  seem 

to  indicate  that  the  cross-section  for  a  typical  exothermic  charge  transfer 

of  this  type  at  about  1  e.  v.  particle  energy  is  of  the  order  of  a  few  times 
2  —  1 6  Z 

(i,  e.  ,  about  3x10  cm  ),  which  is  not  much  smaller  than  the  gas 
kinetic  (momentum  transfer)  cross-section. 

Charge  exchange  reactions  of  the  type  (18-b),  which  involve  a  simultane¬ 
ous  charge  transfer  and  an  explicit  re-arrangement  of  atoms  during  a  col¬ 
lision,  include  the  following  specific  reactions. 


+ 

0 

NO 

+ 

N+ 

- 

2.  2  e,  V. 

(31 -a) 

+ 

0 

=5=5: 

NO"^ 
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N 
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3.  0  e, V. 

(31  -b) 
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NO 
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N+ 

- 

4.  2  e.  V. 

(31-c) 
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NO"^ 
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(31 -d) 
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O2 

=5=t 

NO 
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4,  5  e,  V, 

(32-a) 
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+ _ ^ 

NO 

+ 

no'*' 
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0.  9  e.  V. 

(32-b) 

02'' 
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NO 

+ 
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- 
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(35-a) 

02'' 
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=5=5: 

NO"^ 
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^2 

+  N'*' 

NO  + 

O"*”  +  2.  3  e.  V. 

(35-c) 

^2 

+  N'*' 

NO"''  + 

O  +  6.  7  e.  V. 

(35-d) 

While  no  quantitative  information  appears  to  be  available  in  regard  to  the 

rate  constants  for  these  reactions,  the  similarity  between  these  reactions 

and  the  neutral  atom  exchange  reactions  (4),  (5),  and  (6)  suggests  that  their 

corresponding  exothermic  rate  constants  may  be  of  comparable  magnitude. 

This  would  again  place  the  exothermic  cross-section  for  these  reactions  at 
2 

about  ^  a^. 

The  third  type  of  charge  exchange  reaction  (18-c),  which  involves  a 
simultaneous  charge  exchange  with,  and  ionization  of,  the  neutral  par¬ 
ticle,  would  be  energetically  possible  only  when  the  sum  of  the  first  and 
second  ionization  potentials  of  Y  is  comparable  to  the  first  ionization  po¬ 
tential  of  X.  This  situation  is  not  expected  to  occur  in  the  dissociated  air 
mixture.  * 

D.  Electron  Attachment 

Atomic  and  molecular  oxygen  are  known  to  have  considerable  electron 
affinity.  ^9-41  Thus,  the  attachment  processes. 
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which  are  typical  examples  of  reactions  (19-a),  (19-b),  and  (19-c),  often 
become  the  dominant  rate-determining  processes  for  the  removal  of  free 
electrons  in  many  ionospheric  and  gas-discharge  situations.  ^9,  38,42 
However,  in  view  of  the  fact  that  the  binding  energies  for  and  for  O" 

=i=Note  that  even  for  alkali-earth  elements,  such  as  Ba  and  Ca,  which  may 
be  present  in  small  amounts  as  impurity  in  the  gas  sample,  the  sum  of  the 
first  and  second  ionization  potentials  is  still  greater  than  15  e.  v. 
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are  only  about  0.  5  and  1.  5  e.  v.  ,  respectively,  the  concentrations  of  O2 

and  O”  ions  would  always  be  negligibly  small  in  comparison  with  the  elec- 

o 

tron  concentration  whenever  the  gas  temperature  exceeds  about  1,500  K. 

Even  though  the  effect  of  these  attachment  processes  on  the  ionization 
history  is  negligible,  reactions  (42)  and  (43)  may  contribute  significantly 
to  the  dissociation  rate  of  ©2  and  NO  at  high  electron  temperature  and 
density. 

E.  Ionizing  Atomic  and  Molecular  Collisions 

From  general  considerations  of  inelastic  collisions  between  two  atomic 
systems  under  gas-kinetic  conditions,  one  may  conclude  that  ionizing  colli¬ 
sions  of  the  type  represented  by  Eq.  (20)  would  not  be  important  as  long  as 

3 

the  mean  kinetic  energy  per  particle  kT  is  small  compared  with  the 

ionization  potential  of  the  atom  (or  molecule)  being  ionized.  This 

conclusion  can  be  arrived  at  either  from  a  simple  test  of  the  near -adiabatic 
condition,  ** 


- °  >  >  1  (44) 

h  (kT/m^)^/^ 

or  from  the  usual  argument  that  energetic  collisions  between  two  particles 
in  which  the  energy  of  relative  motion  exceeds  the  ionization  potential  are 
relatively  rare. 

While  the  above  conclusion  is  certainly  correct  for  noble  gases  in  the 

2 

range  of  gas  temperature  under  consideration,  we  must  now  consider  the 
possibilities  that  the  chemical  energy  in  the  colliding  atomic  systems  may 
come  into  play  so  that  the  required  amount  of  energy  transfer  from  relative 
translational  motion  to  electronic  excitation  may  be  considerably  less  than 
the  ionization  potential  of  the  molecular  complex  formed  by  the  collision. 


*Even  though  collisional  detachment  by  neutral  atoms  and  molecules  (see 
Ref.  41)  may  not  be  sufficiently  rapid  to  keep  the  negative  ions  in  quasi - 
equilibrium  with  the  free  electrons  and  the  neutral  particles,  there  is  no 
reason  why  collisional  detachment  by  electrons  should  not  be  sufficiently 
rapid  to  maintain  such  a  quasi- equilibrium  condition  behind  the  shock  front. 

**See  for  example,  Massey  and  Burhop,  loc.  cit.  ,  Chapt.  II. 
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This  may  lead  to  a  violation  of  the  near -adiabatic  condition  (44)  for  a  con¬ 
siderable  fraction  of  all  the  molecular  encounters  even  at  gas  temperatures 

4  o  .43 

of  the  order  of  10  K.  Thus,  in  the  presence  of  crossing,  or  near¬ 
crossing,  of  the  potential  energy  curves  for  the  colliding  atomic  system 
(or  potential  energy  surfaces  in  the  case  of  collision  between  molecules), 
ionizing  reactions  of  the  types  represented  by  (21)  through  (23-b)  coixld  be 
extremely  rapid. 

For  the  oxygen-nitrogen  system  under  consideration,  we  may  now  make 
a  representative  (but  not  exhaustive)  list  of  specific  reactions  which  belong 
in  this  category  according  to  an  ascending  order  of  their  respective  energy 


of  reaction  AE°; 


N 

+ 

o 

+ 

2.  8  e.  V. 

N 

+ 

N 

+ 

5.  8  e.  V. 

N 

+ 

O2 

+ 

6.  5  e.  V. 

O 

+ 

0 

+ 

6,  9  e.  V. 

O 

+ 

NO 

4- 

7.  9  e.  V. 

N 

+ 

NO 

+ 

7.  9  e.  V. 

X 

+ 

NO 

+ 

9.  3  e.  V. 

O 

+ 

^2 

+ 

11.2  e.  V. 

^2 

+ 

O2 

+ 

11.2  e.  V. 

O 

+ 

O2 

+ 

11.7  e.  V. 

X 

t 

O2 

+ 

12.  1  e.  V. 

X 

+ 

0 

+ 

13.  6  e.  V. 

X 

+ 

N 

+ 

14.  6  e.  V. 

X 

+ 

^2 

+ 

15.  6  e.  V. 

no"*"  + 

e 

(45) 

n/  + 

e 

(46) 

no/  + 

e 

(47) 

0/  + 

e 

(48) 

no/  + 

e 

(49) 

N20+  + 

e 

(50) 

X  + 

NO"^ 

+ 

e 

(51) 

4* 

N2O+  + 

0 

(52) 

NO  + 

no'*’ 

+ 

e 

(53) 

0/  + 

e 

(54) 

X  + 

°2^ 

+ 

e 

(55) 

— 

X  + 

0+ 

+ 

e 

(56) 

X  + 

N+ 

+ 

e 

(57) 

X  + 

^2^ 

+ 

e 

(58) 
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From  an  earlier  study  of  the  ionization  rate  behind  shock  waves  in  diluted 
O2  “  ^2  mixtures,  Lin'^'^  has  made  a  crude  determination  of  the  forward 
rate  constant  for  reaction  (45).  The  velocity-averaged  cross-section  asso¬ 
ciated  with  this  process  as  deduced  from  these  earlier  experiments  is  re¬ 
produced  here  in  Fig.  5.  It  may  be  mentioned  that  a  correction  factor  of 
about  5  has  been  applied  to  the  lower  bound  estimate  presented  in  Ref.  44 
to  account  for  the  effect  of  finite  dissociation  rate,  so  that  these  results 
now  represent  more  closely  the  true  cross-section  for  this  process  within 
the  temperature  range  covered  by  the  experiments.  To  extrapolate  these 

results  to  other  temperatures,  we  shall  proceed  in  the  following  manner: 

45 

According  to  Bates  and  Massey,  the  cross-section  associated  with 
each  energetically  allowed  ionizing  collision  (i.  e.  ,  impact  energy  E  greater 
than  the  activation  energy  E^)  due  to  potential  energy  curve  crossing  is 
expressible  in  the  form 


‘E  >  E 


(59) 


where  g  is  the  ratio  between  the  statistical  weight  associated  with  the  initial 

potential  energy  curve  which  leads  to  the  crossing  point  and  the  sum  of  the 

statistical  weights  associated  with  all  possible  initial  potential  energy  curves 

of  the  colliding  atomic  system;  is  the  distance  between  the  two  centers 

of  mass  of  the  colliding  system  where  the  crossing  occurs;  /n  =  "*■ 

is  the  reduced  molecular  weight  of  the  colliding  system;  E  is  the  impact  energy 

(measured  in  electron-volts);  and  T  .  ,  is  the  auto -ionization  life  time  meas- 

ured  in  seconds)  of  the  molecular  complex  formed  after  the  potential  energy 

curve  crossing,  which  is,  of  course,  the  main  unknown  in  the  problem.  Ac- 

46  47 

cording  to  the  worKs  of  Bransden  and  Dalgarno,  Wu,  and  others,  the 


*The  lower  bound  estimate  presented  in  Ref.  44  was  based  on  the  experiment¬ 
ally  observed  ionization  rate  behind  the  shock  but  assuming  that  the 
mixture  was  everywhere  in  dissociation  equilibrium.  More  recent  finite 
rate  calculations  indicated  that  the  dissociation  equilibrium  assumption  has 
overestimated  the  local  atom  concentrations  by  a  considerable  factor. 
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auto-ionization  time  for  some  typical  doubly- excited  states  of  helium  appears 

-15  -11 

to  lie  in  the  range  10  <  j  <  10  sec.  ,  but  no  quantitative  informa¬ 

tion  appears  to  exist  for  more  complex  atomic  systems. 

The  potential  energy  curves  for  the  lower  electronic  states  of  the  N  +  O 
and  the  N  +  O  systems  are  illustrated  in  Fig-  6.  It  is  seen  that  there  exist 
several  possibilities  for  crossing  of  the  N  +  O  curves  near  the  potential  en¬ 
ergy  minimum  of  the  ^  state  of  NO"*".  One  may  therefore  interpret  the 
experimental  result  shown  in  Fig.  5  as  due  to  such  curve-crossing.  Thus, 

for  the  system  under  consideration,  we  have,  p  =  7.5;R  a  =10"  cm; 

-14^  ° 

and  =  3  e.  v.  If  j  were  much  shorter  than  10  sec,  Eq.  (59)  be¬ 
comes  approximately,  * 

Qe  >  E  ^3x10'^'^  cm^  (59-a) 

Behind  the  translational/rotational  shock  front,  the  velocity  distribution  of 
the  atoms  and  molecules  will  be  approximately  Maxwellian  at  the  local  trans¬ 
lational  temperature  T.  The  velocity-averaged  ionization  cross-section  would 

u  48 
then  be, 

Q  =  3  X  10"^”^  (E^kT)  exp  (-E^/kT)  cm^  (60 -a) 

- 1  4 

On  the  other  hand,  if  T.  j  were  much  longer  than  10  sec.  ,  Eq.  (59)  be- 

A.0  X. 

comes  approximately 

^E  >E  =  6  X  10-^1  (E^/2  j  cm^  (59-b) 

X 

Assuming  j  to  be  independent  of  impact  energy,  the  corresponding 
velocity-averaged  cross-section  becomes, 

■Q  =  3  X  10-^1  (E^^/2  j 


*For  the  N  (4^S)  +  0(^P)  system,  g  would  be  either  l/l8  or  2/9,  depending 
on  whether  the  ^2'*’  the  4 jj- state  is  involved  in  the  crossing.  Therefore, 
a  reasonable  number  to  be  assigned  to  g  would  be  0.  1  when  the  state  in¬ 
volved  is  not  known. 
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If  one  tries  to  fit  the  center  of  the  experimental  points  shown  in  Fig.  5  to 
Eq.  (60-a),  which  has  only  one  adjustable  parameter  (assuming  that  only 
one  crossing  point  was  involved  in  the  experimental  situation),  one  obtains 
E  =  4.  8  e.  V.  ,  and 

X  ’ 

U  =  1. 7  X  10"^'^  T"^  exp  (-55,  700/T)  cm^  (6l-a) 

which  gives  the  steepest  possible  temperature  dependence  (dotted  curve  in 
Fig.  5).  On  the  other  hand,  if  we  let  E^  =  2.  8  e.  v.  ,  which  is  the  lowest 
possible  activation  energy  for  reaction  (45),  fitting  the  center  of  the  same 
set  of  experimental  points  to  Eq.  (60-b)  yields  T  .  .  =  6  x  10  sec,  and 

Q  =  10"^'^  T"^  exp  (-32,  500/T)  cm^  (6l-b) 


which  gives  the  weakest  possible  temperature  dependence  (solid  curve  in 
Fig.  5).  Since  the  limited  temperature  range  and  the  considerable  scatter 
of  the  experimental  points  as  shown  in  Fig.  5  do  not  allow  a  better  deter¬ 
mination  of  the  temperature  dependence  within  he  above  two  possible  ex¬ 
tremes,  we  shall  draw  on  some  other  experimental  evidence  for  such  a 
determination. 

The  forward  rate  constants  (i.  e.  ,  the  mean  cross-section  multiplied  by 
the  mean  thermal  speed  based  on  the  reduced  mass  of  the  two  colliding  atoms) 
corresponding  to  (6l-a)  and  (6l-b)  are,  respectively, 

/C45(a)  =9x10"^  exp  (-55,  700/T)  cm^sec  (62-a) 

/C  45(b)  =  5  X  10”^^  exp  (-32,  500/T)  cm^sec  (62-b) 

Within  the  temperature  range  300  ^  T  ^  30,000°K,  the  equilibrium 
constant  for  reaction  (45)  can  be  approximated  by  the  polynomial  (to  an 
accuracy  of  about  +  10%,  assuming  no  error  in  the  energy  of  reaction), 

K..  =  (1.  4  X  10'®  T  +  1.  2  X  10”^^  +  1.  4  x  lO"^^  T^)  exp  (-32,  500/T) 

(63) 
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Dividing  (62-a)  and  (62-b)  by  (63),  one  obtains  the  corresponding  backward 

rate  constants  for  reaction  (45),  which,  at  low  temperatures  (i.  e.  ,  T  << 

4  o 

10  K),  become  approximately, 

'^R45(°a)  =  exp  (-23,  200/T)  cm^sec  (64-a) 

K  R45(b)  =  3  X  10  ^  T  cm^/sec  (64-b) 


Note  that  the  steep  temperature  dependence  (64-a)  would  rule  out  any  experi¬ 
mental  observation  of  dissociative  recombination  for  NO  at  or  near  room 

temperature,  since  the  numerical  value  of  .  would  then  be  of  the 

38  3  rv^o^a^ 

order  of  lO"  cm  /sec.  On  the  other  hand,  the  weak  temperature  depend¬ 
ence  (64-b)  gives  K  =  6  x  lO"^  cm^sec  at  300°K,  which  is  quite 

compatible  with  the  dissociative  recombination  coefficient  of  about  2x10 

3  49 

cm  /sec  deduced  by  Doering  and  Mahan  from  photolysis  of  NO.  At  a 

-  8  3 

temperature  of  2000°K,  the  same  equation  gives  =  3  x  10  cm  /sec, 

which  is  also  compatible  with  Sugden’s  observation*  that  the  dissociative 
recombination  coefficient  for  NO'*'  at  flame  temperatures  appeared  to  be  of 
the  order  of  lO""^  cm^sec.  Therefore,  these  low  temperature  experiments 
seem  to  support  the  weak- temperature  dependence  (62-b)  and  (64-b)  for  the 
rate  constants  of  reaction  (45). 

A  question  may  be  raised  here  in  regard  to  the  validity  of  comparing  the 
forward  and  backward  rates  for  reaction  (45)  in  the  above  manner,  as  the 
electronic  states  of  the  atoms  and  the  detailed  paths  associated  with  the 
observed  rates  are  not  known.  Actually,  in  order  for  the  above  compari¬ 
son  to  be  valid,  a  sufficient  condition  would  be  that  the  rate  governing 
step  was  in  the  potential  energy  curve  crossing,  and  that  both  the  forward 
and  the  backward  processes  passed  through  the  same  crossing  point.  It 
may  further  be  noticed  that  this  also  happened  to  be  a  necessary  condition 


*T.  M.  Sugden,  discussion  following  Doering  and  Mahan* s  paper ,  loc.  cit. 
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in  order  to  explain  Sugden' s  flame  temperature  result. 

For  reactions  (46)  and  (48),  while  no  quantitative  measurement  of  the 
forward  rate  constants  has  yet  been  made,  the  recombination  rates  for 
electrons  and  N,^,  ions  have  been  a  subject  of  study  for  some  time 

in  gas  discharge.  The  recent  work  of  Kasner,  Rogers,  and  Biondi, 
which  confirmed  the  existence  of  dissociative  recombination  in  the  after 
glow  decay  of  a  microwave  discharge  through  the  use  of  a  mass  spectrom¬ 
eter,  indicated  that  the  recombination  coefficients  for  ions  and  elec- 

trons  and  for  ions  and  electrons  are  (5.  9  +  1. 0)  x  10  cm  /sec  and 

-7  3  " 

(3.  8  +  1.0)  X  10  cm  /sec,  respectively.  Since  the  electron  tempera¬ 
ture  relaxation  time  was  generally  much  shorter  than  the  time  over  which  - 

the  dissociative  recombination  process  was  observed  under  the  experi- 

50 

mental  conditions  described  by  Kasner,  et  al,  the  electron  temperature 

could  not  have  been  much  higher  than  the  gas  temperature  in  the  decaying 

after  glow.  Thus,  the  magnitude  of  the  observed  rate  constants  suggested 

that  the  dissociative  recombination  of  electrons  and  N2^,  ions  also 

did  not  require  any  noticeable  amount  of  activation  energy.  Again,  if  one 

attributes  these  dissociative  recombination  processes  to  the  crossing  of 

potential  energy  curves,  the  crossing  points  must  also  have  occurred  near 

+  + 

the  corresponding  potential  energy  minima  of  the  e  +  N2  and  the  e  + 

system,  and  the  auto -ionization  times  associated  with  the  corresponding 

- 1  2 

and  complexes  must  also  be  of  the  order  of  1 0  second.  Such 
interpretation  of  the  observed  low  temperature  dissociative  recombination 


^Interpretation  of  Sugden' s  result  presents  no  difficulty  if  the  dissociative 
recombination  process  results  in  ground  state  N  and  O  atoms.  On  the  other 
hand,  if  the  process  were  e  +  NO^  (IS"^)  ^  J-  N(^D)  +  0(^P)  +  0.  4  e.  v.  or, 
e  +  MO+(lv:+)  ,  ^  M(^.c!)  +  0(^D)  +  0.  8  e.  v.  the  recombination  process  could 
not  proceed  very  far  at  flame  temperatures  on  account  of  the  very  low  energy 
of  reaction  (which  favors  the  ionized  state)  unless  de- excitation  of  the  N  (^D) 
or  O  (^D)  states  were  sufficiently  rapid.  This  implies  that  unless  the  rate 
governing  step  were  in  the  potential  energy  curve  crossing  and  not  in  the 
collisional  de-excitation  of  metastable  atoms,  the  recombination  rate  at 
flame  temperatures  should  be  very  different  from  that  at  room  temperature. 
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coefficients  also  leads  to  the  following  temperature-dependent  backward 
rate  constants  for  reactions  (46)  and  (48). 


K 


K 


R46 

=  3  X  10"^ 

T-V2 

3/ 

cm  /  sec 

(65) 

R48 

=  2  X  10"^ 

T-3/2 

3, 

cm  /sec 

(66) 

-  3/2 

It  may  be  noted  that  the  T  '  temperature  dependence  for  the  backward 

4 

rate  constants  is  the  same  as  that  suggested  by  Sayers,  but  that  the  numeri¬ 
cal  coefficients  for  ^j^4^  ^R48  shown  are  smaller  than  his  suggested 

values  by  factors  of  about  5  and  1.  5,  respectively.  To  obtain  the  forward 
rate  constants,  one  may  multiply  (65)  and  (66)  by  the  corresponding  equilib¬ 
rium  constants  for  reactions  (46)  and  (48),  which,  in  the  temperature  range 
300  S  T  ^  30,  000°K,  may  be  approximated  by  the  polynomials  (again  to 
an  accuracy  of  about  +  10%,  assuming  no  error  in  the  energy  of  reaction), 

K46  =  (3x  10"®  T  +4x  10"^^  T^+  10‘^^T^  -3x  10"^°  T^)  exp  (-67.  300/T)  (67) 


K^g  S(i.6x10"®T+1.2x10"^^T^+3.5x10'^^T^)  exp (-80,  IOO/t)  (68) 


These  yield, 

^^46  ^  ^46  '^R46 

'^48  ^48  '^R48 


=  9  X  10"^^  T"^/^  (1  +  1.  3  X  lO""^  T 

+  3.  3x  10‘®T^-10‘^^T^)exp(-67,  300/T)  cm^sec 

(69) 

=  3.  2  X  10'^^  (1  +  7.  5  X  10'^  T 

+  2.  2  X  10'®  T^)  exp  (-80,  IOO/T)  cm^/  sec 

(70) 


For  the  molecule-molecule  collision  processes  (53),  (55)  and  (58),  one 

may  deduce  the  corresponding  forward  rate  constants  from  the  recent  molec- 

51  52  53 

ular  beam  results  of  Utterback  and  Miller  and  Utterback.  ’  The  mono¬ 
chromatic  (i.  e.  ,  single-energy)  ionization  cross-sections  for  N2  -  N2, 

N2  -  ©2,  and  ©2  -  O.^  collisions  as  reported  in  Refs.  51  and  5  3  are  repro¬ 
duced  here  in  Fig.  7.  By  extrapolating  the  linear  portion  of  these  curves 
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toward  (E  -  E.) — ►  0,  and  integrating  according  to  the  Maxwellian  molecular 
^  48 

velocity  distribution,  one  obtains  the  following  expressions  for  the  forward 
rate  constants  for  reactions  (53),  (55)  and  (58), 

/<53  =  «55  =  4.  5  X  (1  +  3  X  10"^T)  exp  (-141,  OOO/T)  cmVsec 

^2  33  Q /?  ^  3 

=  1.1x10  (1  +  4  X  lO’^^T)  exp  (-141.000/T)  cm  /sec 

(72) 

28  7 /2  8 

=  2.8x10"  (1  +  3  X  10"  T)  exp  (-181,  OOO/T)  cm  /sec 

5o 

(73) 


It  is  interesting  to  note  that  at  temperatures  below  30,  OOO'^K,  the  pre¬ 
exponential  factors  for  '^55  '^55^^*  Kgg^2  shown  above  are  some¬ 
what  smaller  than  those  for  <<43  given  in  Eqs.  (62-b), 

(69)  and  (70),  but  at  T  =  30,000OK,  they  become  comparable.  This  suggests 
that  the  ionization  processes  for  the  O2  -  O^,  ~  ^2  ^2  "  ^2 

observed  by  Utterback  and  Miller  in  the  molecular  beam  experiments  may 
also  be  interpreted  as  due  to  crossing  of  the  potential  energy  hyper  surfaces 
for  the  colliding  molecular  systems.  The  steep  energy  dependence  of  the 
observed  ionization  cross-section  (see  Fig.  7),  however,  indicates  that  the 
corresponding  auto -ionization  time  t.  ,  associated  with  these  molecular 
collisions  is  not  independent  of  impact  energy. 

For  reactions  (47),  (49),  (50),  (51),  (52),  (54),  (56),  (57),  and  for  reac¬ 
tions  (55),  (58)  when  the  colliding  particle  X  is  other  than  N^,  there 
seems  to  be  no  direct  experimental  information  available  in  the  literature. 
However,  in  view  of  the  general  trend  which  has  just  been  established  be¬ 
tween  the  rate  constants  for  those  reactions  where  experimental  informa¬ 
tion  is  available,  it  is  not  unreasonable  to  express  the  forward  rate  con¬ 
stants  for  all  these  reactions  in  the  simple  Arrhenius  form, 

=  A.  (T)  exp  (  -AEp/kT)  (74) 


*H.  S.  W.  Massey  and  E.  H.  S.  Burhop,  loc.  cit.  ,  p.  450. 
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and  to  make  (T)  comparable  to  either  the  pre- exponential  factors  for 

'^45(b)'  '^46'  '^48  ‘^53’  '^55^^  and  /c  53^2. 

It  may  be  pointed  out  that  in  a  recent  treatment  of  the  ionization  relaxa- 

17 

tion  problem  in  hypersonic  air  flows.  Hall,  Eschenroeder ,  and  Marrone, 

54 

and  Eschenroeder  have  made  estimates  of  the  backward  rate  constants 
for  reactions  (55)  through  (58)  according  to  Thomson's  three-body  recom¬ 
bination  theory.  They  have  subsequently  arrived  at  rate  constants  that 
are  in  effect  many  orders  of  magnitude  greater  than  what  we  have  just 
proposed  above  for  these  unknown  reactions  in  the  temperature  range  of 
interest  (i.  e.  ,  T  30,  000°K).  In  view  of  the  fact  that  such  large  rate 
constants  are  not  consistent  with  the  crossing  point  model  (assuming  any 
reasonable  value  for  the  auto -ionization  time  ^  ),  and  the  fact  that  the 
classical  Thomson  theory  is  probably  not  applicable  to  electron-ion  recom- 
bination,  these  high  estimates  may  be  tentatively  disregarded. 


IV.  RESULTANT  IONIZATION  HISTORY  BEHIND  THE  SHOCK 

By  numerical  integration  of  the  various  ionization  rate  equations  accord¬ 
ing  to  the  local  temperature  and  neutral  particle  concentrations  already  de¬ 
termined  (see  Section  II),  one  obtains  the  instantaneous  electron  production 
rate  as  well  as  the  electron  and  positive  ion  densities  as  functions  of  distance 
behind  the  shock.  The  specific  ionization  rate  due  to  the  various  groups  of 
processes  is  plotted  in  Fig.  8  for  three  typical  shock  velocities.  The  initial 
air  density  chosen  for  these  numerical  examples  is  Pj^  =  0.  02  mm  Hg  at  room 
temperature,  but  to  show  the  density  dependence,  the  specific  ionization  rate 
at  pj^  =  2  mm  Hg  is  also  plotted  in  Fig.  8(b).  The  specific  ionization  rate  is 
defined  here  as  the  absolute  value  of  the  time  rate  of  change  of  the  normalized 
electron  density  (following  the  mass  motion  of  the  fluid)  multiplied  by  the  time 
it  takes  the  shock  to  traverse  one  upstream  mean  free  path,  (which  is 

also  roughly  the  mean  binary  collision  time  behind  the  shock).  The  curve 
labeled  "atom-atom"  refers  to  the  sum  of  the  contributions  from  reactions 
(45),  (46)  and  (48),  which  are  calculated  according  to  equations  analogous 


*H.  S.  W.  Massey  and  E.  H.  S.  Burhop,  loc.  cit.  ,  p.  623. 
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to  Eq.  (7 -a),  using  backward  rate  constants  given  by  (64 -b),  (65),  (66)  and 
corresponding  forward  rate  constants  given  by 

Kj=Kj/CRj  (75) 

where  Kj  is  the  equilibrium  constant  for  the  particular  reaction  (j).  It  may 
be  noted  that  for  reaction  (45),  the  forward  rate  constant  so  specified  be¬ 
comes  somewhat  greater  than  that  corresponding  to  Eq.  (62 -b)  at  high  tem¬ 
peratures  (by  a  factor  of  about  2.  4  at  10,  000°K;  and  about  5.  6  at  20,  000°K). 
This  method  of  extrapolating  to  higher  temperatures  is  probably  more 

realistic  than  a  direct  application  of  Eq.  (62-b)  in  view  of  the  noted  steep 
temperature  dependence  for  ^53»  ^55  ^  ,  and  .  The  curve  labeled 

"Molecule-Molecule  and  Atom -Molecule"  refers  to  the  sum  of  the  contribu¬ 
tions  from  reaction  (47)  and  reactions  (49)  through  (58).  The  contributions 
from  reactions  (53),  (55)  and  (58)  are  calculated  according  to  the  forward 
rate  constants  given  by  Eqs.  (71),  (72)  and  (7  3);  while  the  contributions  from 
all  other  reactions  in  this  group  are  calculated  according  to  Eq.  (74),  using 
a  pre- exponential  factor  given  by  the  geometrical  mean  of  those  for 
©2  ^2 

46’  48’  53*  ^  55  *  58  ’ 

The  curve  labeled  " Electron  Impact"  is  calculated  according  to  Eq.  (26), 

2 

talcing  iX  ^ '^iX^  ~  ^^’o  species.  The  curve  labeled"  Photo - 

ionization"  is  calculated  according  to  Eq.  (29),  with  a  simple  correction 
factor  exp  [-p  (x-L)  ]  to  allow  for  attenuation  of  the  far  u.  v.  intensity  at 

large  distances  behind  the  shock,  using  an  absorption  coefficient  given  by 
P  =  10"^^Sn„  cm“^. 

In  these  specific  rate  calculations,  the  charge  exchange  reactions  (30) 

through  (39)  and  (31 -a)  through  (35-d)  are  allowed  to  proceed  simultaneously 

with  the  electron  production  processes.  The  rate  constants  for  these  20 

charge  exchange  reactions  are  based  on  exothermic  charge  transfer  cross- 
2 

section  of  3  ttq^  ,  as  explained  in  Section  III-C. 

According  to  these  numerical  results,  the  most  important  group  of  ioni¬ 
zation  processes  appears  to  be  the  atom-atom  reactions.  This  is  to  be 
followed  in  an  approximately  decreasing  order  of  importance  by  photoioniza¬ 
tion,  electron  impact,  atom-molecule,  and  molecule-molecule  collisions, 
depending  on  the  shock  velocity.  The  predominance  of  the  atom-atom 


-30- 


processes  is  most  pronounced  at  low  shock  velocities  on  account. of  their 
advantage  in  activation  energy.  The  electron-impact  processes,  on  the 
other  hand,  make  negligible  contributions  at  shock  velocities  below  5  km/ 
sec.  ,  but  may  become  the  predominant  processes  at  shock  velocities  greater 
than  10  km/sec. 

Comparing  the  results  plotted  in  Fig.  8(b)  for  the  two  values  of  p^^  which 
differ  by  a  factor  of  100,  it  is  seen  that  these  specific  ionization  rates  are 
relatively  density-independent.  This  is  due  to  the  predominance  of  the  two- 
body  processes  in  the  range  of  air  density  under  consideration. 

Among  the  atom-atom  ionization  processes,  the  N-O  reaction  (45)  appears 
to  be  the  major  contributor  within  the  range  of  shock  velocity  under  consider¬ 
ation.  This  is  illustrated  in  Fig.  9,  where  the  partial  contributions  from  the 
various  atom-atom,  atom-molecule,  and  molecule -molecule  processes  are 
plotted  for  the  case  =  6.  9  km/sec  and  pj^  =  0.  02  mm  Hg.  The  predomi¬ 
nance  of  reaction  (45)  over  the  other  reactions  is  more  pronounced  at  lower 

shock  velocities,  but  at  U  =9  km/sec  the  contribution  from  the  N-N  reac- 

s 

tion  (46)  becomes  comparable  to  that  from  the  N-O  reaction.  It  may  be 
noted  that  the  curves  (45-R),  (46-R)  and  (48~R),  as  well  as  the  curves 
labeled  "atom -atom,  (d  [e]/dt  <  0)"  in  Figs.  8(a),  (b)  and  (c),  refer  to 
the  rate  of  disappearance  of  free  electrons  due  to  the  excess  of  the  backward 
rate  (i.  e.  ,  dissociative  recombination)  over  the  corresponding  forward  rate. 

The  normalized  electron  and  ion  densities  are  plotted  in  Figs.  10(a),  (b), 
and  (c)  as  functions  of  the  normalized  distance  for  the  same  shock  conditions 
as  those  illustrated  in  Figs.  8(a),  (b),  and  (c).  It  is  seen  that  the  electron 
density  generally  builds  up  monotonically  until  the  atom-atom  reactions  run 
into  reverse.  Beyond  this  point,  the  electron  density  is  continuously  re¬ 
duced  by  dissociative  recombination  as  the  atom-atom  reactions  try  to 
follow  the  decreasing  translational  temperature  toward  the  final  equilibrium 
level. 

Among  the  positive  ions,  NO^  is  expected  to  be  the  most  abundant  species 
up  to  fairly  high  shock  velocity.  However,  at  about  9  km/sec  velocity,  the 
atomic  ions  n"*”  and  O"*”  may  become  more  abundant  than  NO  at  some  distance 
behind  the  shock.  This  is  to  be  expected  at  high  temperatures  in  the  pres¬ 
ence  of  quasi- statistical  equilibrium  between  the  various  positive  ions 
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brought  about  by  the  rapid  charge  exchange  reactions.  It  may  be  mentioned 
that  in  the  numerical  calculations,  radiative  recombination  between  the  atomic 
ions  and  electrons  has  been  neglected  so  that  the  production  and  destruction 
of  these  ions  were  completely  governed  by  the  charge  exchange  reactions. 

The  effect  of  changing  the  charge  exchange  reaction  rate  on  the  resultant 
electron  and  ion  density  distributions  behind  the  shock  is  illustrated  in  Fig. 

11  for  the  case  =  6.  9  km/sec,  and  p^  =  0.  02  mm  Hg.  In  this  figure,  the 
resultant  distributions  of  [e]  ,  [n"^]  ,  [o"^]  ,  and  ([02^]  obtained  by  increas¬ 
ing  and  decreasing  the  adopted  exothermic  charge  exchange  cross-section 
2  ■ 

=  3  7ra  by  a  factor  of  100  are  plotted  for  comparison  with  those 
obtained  from  the  adopted  cross-section.  The  distributions  of  [NO^J  •  [^2^ 
showed  much  smaller  changes  than  those  for  the  other  positive  ions  and 
hence  are  not  plotted.  It  is  seen  that  the  maximum  effect  on  fe]  ,  which  oc¬ 
curs  in  the  neighborhood  of  the  peak  electron  density,  is  only  of  the  orde^r 
of  10  percent.  For  the  positive  ions;  the  small  effect  of  increasing 
by  a  factor  of  100  indicates  that  the  adopted  cross-section  was  already  suf¬ 
ficiently  large  to  keep  the  various  ions  in  quasi- statistical  equilibrium.  In  fact, 
the  relatively  mild  changes  in  Cn+]  ,  [o+],  and  Co2'*'3  obtained  by  decreasing 
Qch  Qjj  factor  of  100  shows  that  it  would  indeed  be  difficult  to  keep  the 

various  ions  away  from  the  state  of  quasi- statistical  equilibrium  behind  the 
shock. 

Before  leaving  this  section,  we  shall  now  examine  briefly  the  effects  of 

upstream  photoionization  and  electron  diffusion  on  the  resultant  ionization 

profile  behind  the  shock.  From  the  experimental  observation  of  Camm, 

35 

Kivel,  Taylor,  and  Teare,  it  is  quite  clear  that  sufficiently  energetic 

photons  could  escape  from  the  translational/rotational  shock  front  to  photo- 

ionize  molecular  oxygen  upstream.  The  p-oduction  of  free  electrons  ahead 

of  a  plane  shock  wave  by  this  "precursor  radiation"  has  been  treated  in 

some  detail  by  Hammerling  in  Ref.  37.  There  is  also  evidence  of  such 

56 

photoionization  effect  being  observed  during  satellite  re-entry.  In  the 
present  discussion,  however,  we  need  only  to  estimate  the  cumulative  ef¬ 
fect  of  such  upstream  photoionization  on  the  resultant  ionization  history 
behind  the  shock. 
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The  electron  density  distribution  behind  the  shock  with  and  without  con¬ 
sidering  upstream  photoionization  is  shown  in  Fig.  12  for  the  case  U  =  6.  9 
km/sec  and  p^  =  0.  02  mm  Hg.  The  solid  curve,  with  no  upstream  photo¬ 
ionization  and  identical  to  that  shown  in  Fig.  10(b),  was  obtained  by  inte¬ 
grating  the  specific  ionization  rates  behind  the  shock,  assuming  [e]  = 

0  at  X  =  0.  The  dotted  curve,  on  the  other  hand,  was  obtained  by  integrat¬ 
ing  the  same  set  of  rate  equations,  but  letting  j[e]  =  lO”^  at  x  =  0.  This 
initial  value  of  [e]  was  obtained  by  a  balance  between  the  outgoing  ionizing 

O  C 

photon  flux  (hv  S:  12.  1  e.  v.  )  observed  by  Camm,  et  al  [see  Eq.  (28)] 
and  the  flux  of  ion  pairs  convected  back  into  the  translational/rotational 
shock  front  per  unit  area  due  to  the  mass  motion  of  the  fluid.  It  is  seen 
that  even  though  the  neglect  of  upstream  photoionization  would  lead  to  seri¬ 
ous  error  in  the  very  early  part  of  the  ionization  history,  by  the  time  [e] 
reaches  about  10  percent  the  peak  value  behind  the  shock,  the  error  involved 
becomes  smaller  than  a  factor  of  two  in  this  case.  It  may  be  pointed  out, 
however,  that  such  error  would  be  somewhat  more  severe  at  higher  shock 
velocities  on  account  of  the  increasing  importance  of  the  electron  impact 
ionization  process  with  shock  velocity  [compare  Figs.  8(b)  and  (c)]  . 

In  regard  to  the  effect  of  electron  diffusion,  one  may  note  that  even  at 
Pj^  =  0.  02  mm  Hg,  the  Debye  shielding  distance' 

2  .1/2 
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=  (kT/47rn^  e^  ) 


(76) 


corresponding  to  the  electron  density  profiles  illustrated  in  Fig.  10(a) 
through  12  is  always  less  than  5  x  10”^cm.  This  is  sufficient  to  insure 
ambipolar  diffusion  ’  and  charge  neutrality  down  to  a  distance  scale 
smaller  than  0.  2  J  .  Since  ambipolar  diffusion  is  not  much  more  rapid 
than  momentum  (or  viscosity)  diffusion,  it  can  again  be  neglected  outside 
the  translational/rotational  relaxation  zone. 


V.  COMPARISON  WITH  EXPERIMENT 

The  theoretical  results  obtained  above  may  now  be  used  for  comparison 
with  the  experimental  observations  of  Lin,  Fyfe  and  Neal  reported  in  the 
complementary  paper  (Ref.  1).  Since  no  direct  experimental  determination 
of  the  positive  ion  distribution  has  yet  been  made,  we  shall  restrict  our 
_  attention  to  the  electron  density  distribution. 
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The  theoretical  electron  density  profile  in  physical  scales  for  the  case 
U  =  6.  9  km/sec,  p.  =  0.  02  lYim  Hg  is  plotted  in  Fig.  13  for  comparison 

S  X 

with  the  electron  density  profile  obtained  from  four  separate  experiments 
taken  at  =  6.  9  +  0.  05  km/sec,  at  the  same  initial  air  pressure.  Curve 

(1)  is  identical  to  the  normalized  curve  fe]  shown  in  Fig.  10(b).  Curves 

(2)  and  (3)  are  reference  curves  obtained  by  changing  the  absolute  rate 

constant  for  the  predominant  ionizing  reaction  (45)  by  a  factor  of  three  in 

-  3  -  3/2  -  3 

either  directions  (i.  e.  ,  taking  =  9x10  T  and  =  1x10 

respectively),  while  holding  all  other  rate  constants  the  same  as 
those  specified  in  the  preceding  section.  The  effect  of  upstream  photoioniza¬ 
tion  is  hardly  noticeable  in  this  linear  plot  since  the  initial  electron  density 
amounts  to  only  about  7  x  10^^  electrons/cm  just  ahead  of,  and  4  x  10 
electrons/cm^  just  behind  the  translational/rotational  shock  front. 

It  is  seen  that  the  front  part  of  the  experimentally  observed  profile  ap¬ 
pears  to  fall  between  curve  (1)  and  curve  (2)  but  the  rear  part  of  the  observed 
profile  seems  to  fall  systematically  above  the  theoretical  curves.  Even  though 
the  apparent  discrepancy  of  the  order  of  20  percent  is  still  within  the  probable 
inaccuracy  of  the  experiment  in  absolute  electron  density  determination,  the 
distinct  difference  in  slope  between  the  theoretical  and  experimental  profiles 
may  be  significant.  However,  in  view  of  the  fact  that  no  reasonable  adjust¬ 
ment  of  either  the  ionization  or  the  chemical  rate  constants  could  reduce 
the  theoretical  rate  of  decrease  of  electron  density  behind  the  peak  to  the 
experimentally  observed  rate,  it  is  believed  that  the  slower  observed  rate 

of  decrease  of  D  at  large  distances  behind  the  shock  front  was  due  to  the 
®  6 
slight  deceleration  of  the  shock  wave  in  the  experimental  situation.  In  a 

decelerating  shock  wave,  the  fluid  elements  further  away  from  the  shock 

front  are  those  that  have  been  heated  by  a  stronger  shock  wave,  so  that  the 

resultant  perturbation  is  in  the  direction  of  reducing  the  rate  of  decrease 

of  translational  temperature  behind  the  shock. 

The  theoretical  peak  electron  density  and  the  maximum  electron  density 

gradient  behind  the  shock  at  initial  air  pressure  pj^  =  0.  02  mm  Hg  are 

plotted  in  Figs.  14  and  15,  as  functions  of  the  shock  velocity  (or  the  shock 

Mach  number  M  =  U  /a.,  where  a.  is  the  speed  of  sound  for  air  at  room 

S  S  1  1 

temperature).  This  may  be  compared  with  the  corresponding  measured 
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quantities  presented  in  Figs.  8  and  10  of  Ref.  1.  The  theoretical  curves  (1), 
(2)  and  (3)  are  again  based  on  the  electron  density  profiles  calculated  accord¬ 
ing  to  /C  =  3  X  10"^  9  X  10"^  and  1  x  10‘^  respec¬ 

tively.  The  electron  density  corresponding  to  the  final  equilibrium  value  far 
behind  the  shock  is  also  plotted  for  comparison,  it  is  seen  that  the  amount 
of  " ionization  overshoot"  observed  in  the  experiment^’  ^  was  perfectly  reason¬ 
able  according  to  the  present  theoretical  model,  even  though  the  relative  in¬ 
sensitivity  of  the  peak  electron  density  to  the  ionization  rate  mcikes  this  quan¬ 
tity  a  poor  measure  of  the  rate  constant.  The  maximum  electron  density 
gradient  shown  in  Fig.  15,  on  the  other  hand,  is  sensitive  to  the  ionization 
rate  and  shows  a  clear  preference  for  curves  (1)  and  (2)  over  curve  (3).  It 
is  also  interesting  to  note  that  according  to  these  calculations,  the  "ioniza¬ 
tion  overshoot"  phenomenon  is  expected  to  disappear  at  shock  velocities  be¬ 
low  4  km/sec  and  above  9  km/sec. 

The  theoretical  "ionization  rise  distance"  and  " ionization  incubation  dis¬ 
tance"  are  plotted  in  Figs.  l6  and  17  for  comparison  with  the  corresponding 
experimental  results  presented  in  Figs.  11  and  12  of  Ref.  1.  These  quan¬ 
tities  are  defined  here  as  the  peak  electron  density  divided  by  the  maximum 
electron  density  gradient;  and  the  distance  for  the  local  electron  density  to 
reach  one -half  the  peaik  value  behind  the  shock,  respectively.  The  curves 
(1),  (2)  and  (3)  have  the  same  meaning  as  before.  It  may  be  noted  that  the 
experimental  points  presented  include  measurements  both  from  the  micro- 
wave  reflection  probe  (hollow  points)  and  from  the  magnetic  induction  probe 
(solid  points).  The  microwave  probe  measured  the  electron  density  directly, 
but  the  magnetic  probe  measured  only  the  d.  c.  electrical  conductivity,  which 
is  related  to  the  electron  density  by 


CT 
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(77) 


However,  within  the  range  of  shock  strength  covered  by  the  experiments, 

(i.  e.  ,  U  ^7  km/sec),  the  electron  collision  frequency  v  calculated  accord- 
ing  to  the  theoretical  temperature,  density  and  chemical  composition  pro¬ 
files  (Section  II)  and  the  published  values  for  the  electron  momentum  trans- 

58 

fer  cross-sections  for  the  various  species  does  not  vary  by  more  than  a 
factor  of  2  with  distance  behind  the  shock  front,  while  the  electron  density 
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changes  by  orders  of  magnitude.  Therefore,  the  measurements  obtained  by 
both  methods  may  be  compared  directly  with  the  theoretical  curves.  It  is 
seen  that  the  experimental  points  are  in  reasonably  good  agreement  with  the 
geometrical  mean  of  curves  (1)  and  (2). 

The  initial  pressure  dependence  of  the  peak  electron  density  and  of  the 
maximum  electron  density  gradient  behind  the  shock  reported  in  Ref.  1  may 
be  interpreted  as  a  direct  consequence  of  the  predominance  of  binary  colli¬ 
sions  in  the  forward  rate  processes  within  the  range  of  gas  density  under 
consideration.  This  is  perfectly  consistent  with  the  theoretical  results  [see 
Figs.  1,  2,  8(b)  and  10(b)] . 

VI.  DISCUSSIONS 


From  the  foregoing  results,  it  is  seen  that  all  the  experimental  observa¬ 
tions  reported  in  Ref.  1  can  be  satisfactorily  interpreted  according  to  the 
present  theoretical  model.  Even  though  many  of  the  ionization  rate  constants 
are  still  uncertain,  comparison  of  the  specific  rates  based  on  tentative  rate 
constants  [see  Figs.  8(a),  (b),  (c)  and  ll]  showed  that  only  a  few  of  them 
are  pertinent.  Therefore,  the  number  of  rate  constants  that  can  be  adjusted 
to  fit  the  experimental  results  is  actually  quite  limited.  In  fact,  the  appar¬ 
ent  predominance  of  the  atom-atom  reactions  throughout  the  range  of  shock 
velocity  under  consideration  showed  that  only  through  an  adjustment  of  the 
rate  constant  for  reaction  (45),  and  to  a  lesser  extent,  the  rate  constants  for 
reactions  (46)  and  (48),  can  one  effectively  alter  the  ionization  profile  behind 
the  shock. 


The  reasonably  close  (i.  e.  ,  within  a  factor  of  2)  agreement  between  the 
experimental  points  and  the  theoretical  curve  (1)  shown  in  Figs.  13  through 
17  suggests  that  the  assignment  of  K  and  K  according  to  Eqs. 

(64-b),  (65)  and  (66)  is  already  quite  satisfactory.  Even  though  an  upward 
revision  of  K  by  a  factor  of  2  would  bring  the  theoretical  curve  to  near 
perfect  agreement  with  the  experimental  points,  such  a  fine  adjustment  may 
not  be  justified  at  this  time  in  view  of  the  existence  of  other  uncertainties 


in  the  theoretical  model  which  could  be  of  comparable  magnitude.  These 
include,  especially,  uncertainties  in  the  chemical  composition  profile 
caused  by  uncertainties  in  the  high  temperature  dissociation  rate  constants,  ’ 
and  in  the  model  of  coupling  between  vibrational  relaxation  and  dissociation. 

The  theoretical  results  presented  here  may  be  extended  over  a  much  wider 
range  of  initial  air  density  without  difficulty.  The  upper  density  limit  is 
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reached  when  the  rates  of  dissociation  and  ionization  due  to  three -body 
collisions  (comparing. only  the  forward  rates)  become  comparable  to  those 
due  to  binary  collisions.  This  is  not  expected  to  occur  at  densities  much 
below  liquid  air  density.  For  extension  to  lower  densities,  while  there 
seems  to  be  no  theoretical  limit,  a  practical  limit  should  be  set  by  the 
realizability  of  air  samples  without  excessive  contamination  from  the 
vapor  pressure  of  the  container,  or  from  the  environment,  if  the  theoreti¬ 
cal  results  were  to  be  compared  with  experimental  observations. 

The  theoretical  calculation  can  readily  be  extended  to  other  initial  O2/N2 
mixture  ratios,  and  also  to  predissociated  oxygen-nitrogen  mixtures  for 
upper  atmosphere  application.  With  minor  modification,  the  effects  of 
adding  small  amounts  of  argon,  and  other  gases  to  the  mixture 

may  also  be  estimated,  using  rate  constants  of  the  form  given  by  Eq.  (74). 
Contributions  to  the  electron  impact  ionization  rates  by  these  extra  gases 
can  readily  be  included  in  Eq.  (26). 

Extension  of  the  present  theoretical  model  to  a  wider  range  of  shock 
strength,  on  the  other  hand,  represents  a  more  complicated  problem.  As 
already  mentioned  earlier,  extension  of  the  shock  velocity  above  10  km/sec 
not  only  requires  better  knowledge  of  the  electron  impact  ionization  rates 
(Section  III-A),  but  also  necessitates  consideration  of  the  two-way  coupling 
between  the  dissociation  processes  and  the  ionization  processes.  Extension 
of  the  shock  velocity  below  4  km/sec  is  a  lesser  problem,  but  one  must  then 
consider  the  increasing  role  of  impurities  having  low  ionization  potentials 
with  decreasing  gas  temperature. 

One  remark  may  also  be  made  concerning  the  applicability  of  the  present 
one -dimensional  results  to  three-dimensional  problems.  As  is  well  known 
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in  fluid  dynamics,  there  is  a  one-to-one  correspondence  between  the  flow 
across  an  oblique  shock  wave  and  the  flow  across  a  normal  shock  wave. 
Therefore,  the  present  results  are  directly  applicable  to  plane  oblique  shock 
waves  provided  that  only  the  normal  component  of  the  up-stream  flow  velo¬ 
city  is  used  to  determine  the  shock  strength,  and  that  the  distance  scale  re¬ 
fers  to  the  normal  distance  from  the  shock  surface.  With  slight  modifications. 


the  present  model  can  also  be  used  to  follow  the  ionization  history  in  the  three- 
dimensional  inviscid  flow  field  of  hypersonic  objects^^’  and  in  high-enthalpy 
nozzle  flows. 
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Us=  5.0  KM /SEC 
p,  =  0.02  MM  HG 
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Fig.  1  (a)  Typical  temperature  and  density  distribution  behind  a  normal 

shock  wave  in  air  calculated  from  coupled  vibrational  excitation 
and  chemical  reactions.  X/i.  i  is  the  distance  behind  the  trans¬ 
lational/rotational  shock  front  divided  by  the  upstream  mean  free 
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Fig.  1  (b)  Typical  temperature  and  density  distribution  behind  a  normal 

shock,  wave  in  air  calculated  from  coupled  vibrational  excitation 
and  chemical  reactions.  X/JL  is  the  distance  behind  the  trans- 
lational/4'otational  shock  front  divided  by  the  upstream  mean  free 
na  th . 


Fig.  1  (c)  Typical  temperature  and  density  distribution  behind  a  normal 

shock  wave  in  air  calculated  from  coupled  vibrational  excitation 
and  chemical  reactions.  X/JL  is  the  distance  behind  the  trans¬ 
lational/rotational  shock  front  divided  by  the  upstream  mean  free 


Fig.  Z  (a)  Typical  neutral  particle  distribution  behind  a  normal  shock  wave 
in  air  calculated  from  coupled  vibrational  excitation  and  chemical 
reactions.  [X]  is  the  total  number  of  type  X  particles  at  all  in¬ 
ternal  energy  states  per  unit  volume  divided  by  the  equivalent 
number  of  air  molecules  within  the  same  volume. 


-46- 


9014 


UJ 


AIR  SHOCK 

Ugs  6.9  KM /SEC 

-  p,  »  0.02  MM  HG 

- p  e  2  MM  HG 


i  1 1  i-1  _ 

10 


Fig.  2  (b)  Typical  neutral  particle  distribution  behind  a  normal  shock  wave 
in  air  calculated  from  coupled  vibrational  excitation  and  chemical 
reactions.  [X]  is  the  total  number  of  type  X  particles  at  all  in¬ 
ternal  energy  states  per  unit  volume  divided  by  the  equivalent 
number  of  air  molecules  within  the  same  volume. 


Fig.  2  (c)  Typical  neutral  particle  distribution  behind  a  normal  shock  wave 
in  air  calculated  from  coupled  vibrational  excitation  and  chemical 
reactions.  [X]  is  the  total  number  of  type  X  particles  at  all  in¬ 
ternal  energy  states  per  unit  volume  divided  by  the  equivalent 
number  of  air  molecules  within  the  same  volume. 


-48- 


9016 


10^/ T 


Fig.  5  Velocity-averaged  ionization  cross  -  section  for  N-O  collisions 

deduced  from  observation  of  the  ionization  rate  behind  shock  waves 
in  gas  mixtures  containing  a  small  percentage  of  O2  in  N2  (Ref.  44). 
Note  that  a  correction  factor  of  about  5  has  been  applied  to  the  pre¬ 
vious  lower  bound  estimate  to  allow  for  the  finite  dissociation  rates 
of  O2  and  N2  (see  footnote  on  page  24). 
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POTENTIAL  ENERGY,  e.V. 


Potential  energy  curves  for  the  lower  electronic  states  of  the  NO 
and  the  NO"*"  systems.  (  Courtesy  F.  R.  Gilmore,  Rand  Corp.  ) 


(E-Ej),  e.v.,  CENTER  OF  MASS  SYSTEM 


Ionization  cross-sections  for  N2-N2,  N2-O2  and  O2-O2  collisions 
obtained  from  molecular  beam  experiments  by  Utterback  and  Mille 
(  Ref.  51)  and  by  Utterback  (Refs.  52  and  53). 


o 


Fig.  8  (a)  Specific  ionization  rates  for  the  various  groups  of  ionization 
processes  as  functions  of  distance  behind  the  shock  at 
Ug  =  5,0  km/ sec  and  pj  =  0.  02  mm  Hg. 
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Fig.  9  Detail  breakdown  of  specific  ionization  rates  for  the  various  atom- 
atom,  atom-molecule  and  molecule-molecule  reactions  as  functions 
of  distance  behind  the  shock  for  the  case  Ug  =  6.  9  km/ sec  and 
Pj^  =  0.  02  mm  Hg.  The  curves  labelled  (45-R),  (46-R)  and  (48-R) 
indicate  negative  electron  production  rate  due  to  the  excess  of  the 
backward  rate  over  the  forward  rate  for  reactions  (45),  (46)  and  (48) 
respectively. 
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Fig.  10  (b)  Distribution  of  free  electrons  and  various  positive  ions 

behind  the  shock  at  Ug  =  6.  9  km/sec  and  at  two  different 
values  of  the  initial  pressure  p^.  These  results  are  ob¬ 
tained  by  integration  of  the  resultant  net  production  rate 
of  the  corresponding  species  due  to  the  various  ionization 
processes.  The  effects  of  diffusion  and  upstream  photo¬ 
ionization  have  been  neglected,  but  charge  exchange  reactions 
are  allowed  to  proceed  simultaneously  with  the  production 
processes,  assuming  a  universal  exothermic  charge  exchange 
cross-section  of  Sir  a?  . 


Distribution  of  free  electrons  and  various  positive  ions 
behind  the  shock  at  Ug  =9.0  km/ sec  and  =  0.  02  mm  Hg 
These  results  are  obtained  by  integration  of  the  resultant 
net  production  rate  of  the  corresponding  species  due  to  the 
various  ionization  processes .  The  effects  of  diffusion  and 
upstream  photoionization  have  been  neglected,  but  charge 
exchange  reactions  are  allowed  to  proceed  simultaneously 
with  the  production  processes,  assuming  a  universal  exo¬ 
thermic  charge  exchange  cross-section  of  3  tt  . 
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Fig,  11  Effects  of  varying  the  universal  exothermic  charge  exchange  cross- 
section  by  a  factor  of  100  in  either  direction  on  the  distributions  of 
free  electrons  and  N+,  0+  and  O2+  ions  behind  the  shock  for  the  cas 

®  Pi  I  fli-  The  corresponding  effects 

on  the  distribution  of  NO  and  N2+  ions  not  shown  are  between  those 
01  the  electrons  and  of  the  other  three  positive  ions. 
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Fig.  14  Comparison  between  calculated  and  observed  peak  electron  density 
behind  the  shock  as  a  function  of  shock  velocity  at  pj^  =  0.  OZ  mm  Hg. 
The  theoretical  curves  (1),  {Z)  and  (3)  have  the  same  meaning  as  in 
Fig.  13.  The  curve  (ne)equil.  gives  the  theoretical  electron  density 
corresponding  to  the  final  adiabatic  equilibrium  state  far  behind  the 
shock.  The  experimental  points  are  taken  directly  from  Fig.  8 
of  Ref.  1. 


.  ELECTRON  DENSITY  GRADIENT, 

'dx 


U3,  KM/SEC 


Fig.  15  Comparison  between  calculated  and  observed  maximum  electron 
density  gradient  behind  the  shock  as  a  function  of  shock  velocity 
at  Pi  =  0.  02  mm  Hg.  The  theoretical  curves  (1),  (2)  and  (3)  hav« 
the  same  meaning  as  in  Fig.  13.  The  experimental  points  are  ta 
directly  from  Fig.  10  of  Ref.  1. 
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Fig.  17  Comparison  between  calculated  and  observed  ionization  incubation 
distance  behind  the  shock  as  a  function  of  shock  velocity  at  = 
0.02  mm  Hg.  The  theoretical  curves  (1),  (2)  and  (3)  have  the  same 
meaning  as  in  Fig.  13.  The  experimental  points  are  taken  directly 
from  Fig.  12  of  Ref.  1. 
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